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Rare A^ -^ Ai^i decay is investigated in framework of general two Higgs doublet 
model, in which a new source of CP violation exists (model III). The polarization 
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Oh. asymmetry is a very sensitive tool for establishing model III. 
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1 Introduction 

Rare decays, induced by flavor-changing neutral current (FCNC) b -^ s{d) transitions, 
provide testing grounds for the standard model (SM) at loop level and can give valuable 
information about the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vtd, Vts, Vtb, 
etc. In addition, the study of rare decays can pave the way for establishing new physics 
beyond SM, such as two Higgs doublet model (2HDM), minimal supersymmetric extension 
of the SM (MSSM). Most important of all, study of the h — >• s{d)i~^i~ decays is expected to 
be one of the most reliable quantitative tests of FCNC. This transition has been extensively 



investigated in framework of the SM, 2HDM and MSSM ll|]-[||]. 

The matrix element of the b —>■ si~^£^ contains terms describing the virtual effects 
induced by ti, cc and uu loops which are proportional to VtbV^*, VcbV*^ and VubV*^, re- 
spectively. Using unitarity of the CKM matrix and neglecting Vub V*^ in comparison to 
Vtb Vfl and Vcb V*g, it is obvious that the matrix element for the b -^ si~^i~ involves only one 
independent CKM factor Vtb V^* so that CP-violation in this channel is strongly suppressed 
in the SM. 

The present work is to devoted studying the exclusive A^ -^ Ai'^i^ decay, which at 
quark level is described by 6 — *> s£^£~ transition, in context of the general two Higgs 
doublet model in which a new source for CP violation exists. 

2HDM model is one of the simplest extension of the SM, which contains two complex 
Higgs doublets, while the SM contains only one. In general, in 2HDM the flavor changing 
neutral currents (FCNC) that appear at tree level, are avoided by imposing an ad hoc 
discrete symmetry [|T7|. One possible solution to avoid these unwanted FCNC at tree level 
is that all fermions couple to only one of the above-mentioned Higgs doublets (model I). 
The other possibility is the coupling of the up and down quarks to the flrst and second Higgs 
doublets, with the vacuum expectation values V2 and Vi, respectively (model II). Model II is 
more attractive since its Higgs sector coincides with the ones in the supersymmetric model. 
The strength of couplings of fermions with Higgs flelds depends on tan/3 = f2/fi, which 
is the free parameter of the model. The new experimental results of CLEO and ALEPH 
Collaborations [^8], |19| on the branching ratio b ^ s'y decay impose strict restrictions on the 
charged Higgs boson mass and tan/5. Recently, the lower bound on these parameters were 



determined from the analysis of the 6 ^ 57 decay, including NLO QCD corrections [^, ^ 
The phenomenological consequence of a more general model in 2HDM, namely, model HI, 
without discrete symmetry has been investigated in [|22[-||2^. In this model FCNC appears 
naturally at tree level. However, the FCNC's involving the flrst two generations are highly 
suppressed, as is observed in the low energy experiments, and those involving the third 
generation is not as severely suppressed as the flrst two generations, which are restricted 
by the existing experimental results. 

In this work we assume that all tree level FCNC couplings are negligible. However 
even with this assumption, the couplings of fermions to Higgs bosons may have a complex 
phase e"^. In other words, in this model there exists a new source of CP violation that is 
absent in the SM, model I and model II. The effects of such an extra phase in the 6 — > S7 
and b — > di~^i', B — > 7r£+£~ and B -^ pi~^i~ decays were discussed in [^ |2^ and p7| , 
respectively. The constraints on the phase angle in the product Xu^bb of Higgs-fermion 
coupling (see below) imposed by the neutron electric dipole moment, B^ — B^ mixing, po 



parameter and Rb is discussed in 

The paper is organized as follows: In Section 2 we present the necessary theoretical 
framework and the branching ratio, CP-violating effects in the partial widths for the above- 
mentioned exclusive decay channels are studied. Section 3 is devoted to the numerical 
analysis and concluding remarks. 



2 Theoretical calculations for the A5 -^ Ai^i decay 

Before presenting the theoretical results for A^ — > Ai~^i~ decay, let us remember the main 
essential points of the general Higgs doublet model (model III). Without loss of generality 
we can work in a basis such that only the first doublet generates all the fermion and gauge 
boson masses, whose vacuum expectation values are 
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In this basis the first doublet 0i is the same as in the SM, and all new Higgs bosons result 
from the second doublet 02, which can be written in the following form 




V2 



( V2H+ 



where G~^ and G^ are the Goldstone bosons. The neutral Xi ^iid X2 ^^^ ^'^^ ^^^ physical 
mass eigenstate, but their linear combinations give the neutral H^ and h^ Higgs bosons: 
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H^ cos a — hP sin a 
■- if ° sin a + h^ cos a 



The general Yukawa Lagrangian can be written as 



(1) 



V^'^ and i^ 



jj' , in general, are the non- 



where i, j are the generation indices, = ia2(, 
diagonal coupling matrices, L = (1— 75)72 and R = (l+75)/2 are the left- and right-handed 
projection operators. In Eq. (1) all states are weak states, that can be transformed to the 
mass eigenstates by rotation. In mass eigenstates the Yukawa Lagrangian is 
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where U{T>) represents the mass eigenstates of u, c, t {d, s, b) quarks. In this work, we 
will use a simple ansatz for C,^ '-^ 
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(3) 



assume that Ajj is complex, i.e., Ajj = \Xij\ e*"^. For simplicity we choose ^^' to be diagonal 
to suppress all tree level FCNC couplings, and as a result Xij are also diagonal but remain 
complex. Note that the results for model I and model II can be obtained from model III 
by the following substitutions: 



Xtt = cot P Xbb = — cot P for model I , 

Xtt = cot P Xbb = + tan /3 for model II 



(4) 



and = 0. 

After this brief introduction about the general Higgs doublet model, let us return our 
attention to the b — > si~^i~ decay. The powerful framework into which the perturbative 
QCD corrections to the physical decay amplitude incorporated in a systematic way, is the 
effective Hamiltonian method. In this approach, the heavy degrees of freedom, t quark, 
W^, H^, h^, H^ are integrated out. The procedure is to match the full theory with the 
effective theory at high scale /i = rriw, and then calculate the Wilson coefficients at lower 
fi ~ 0{mb) using the renormalization group equations. In our calculations we choose the 
higher scale as /i = mw, since the charged Higgs boson is heavy enough {mH± > 210 GeV 
20| ) to neglect the evolution from mij± to niw. 



see 



In this work the charged Higgs boson contributions are taken into account and the 
neutral Higgs boson exchange diagram contributions are neglected since Higgs-fermion 
interaction is proportional to the lepton mass. The charged Higgs boson exchange diagrams 
do not produce new operators and the operator basis is the same as the one used for the 
h —>■ si^i~ decay in the SM. Therefore in model HI, the charged Higgs boson contributions 
to leading order change only the value of the Wilson coefficients at m^/ scale, i.e.. 
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The coefficients Cf^^^'^{my/) to the leading order are given by 
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H± 



and sin^ 9^ = 0.23 is the Weinberg angle. It follows from Eqs. (5-8) that among all the 
Wilson coefficients, only Cj involves the new phase angle 0. 
The effective Hamiltonian for the b -^ si'^i' decay is [28-31] 



Grp 
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where Cj are the Wilson coefficients. The explicit form of all operators Oj can be found in 
[28-31]. 

The evolution of the Wilson coefficients from the higher scale /i = mw down to the low 
energy scale /t. = ?77b is described by the renormalization group equation 

where 7 is the anomalous dimension matrix. The coefficient C7 (/i) at the scale /i = ?77b 
in next to leading order (NLO) is calculated in ^^, |21|: 

4:71 

where Cj{mi,) is the leading order (LO) term and C7' (777^) describes the NLO terms, 
whose explicit forms can be found in [^. In our case, the expressions for these coefficients 
can be obtained from the results of ||2l| by making the following replacements: 



\Y\ 



|A«| 



and XY* -^ \\tAbb\e 



i(f> 



In the SM, the QCD corrected Wilson coefficient Cglmh), which enters to the decay amph- 
tude up to the next leading order has been calculated in [28-31]. The Wilson coefficient Cio 
is not modified as we move from /x = mw to fi = nib scale, i.e., Cio^irib) = CfQ^'^{mw)- 
As we have already noted, in model III there does not appear any new operator other than 
those that exist in the SM, therefore it is enough to make the replacement C^*^(mvi/) — *■ 
Cg^^^'^mw) in [28-31], in order to calculate Cg^^^'^ at nib scale. Hence, including the 
NLO QCD corrections, Cg{mb) can be written as: 



Cg{^,) = Cr-'\^,) 



1 + ^^uis) 



+ girhc, s) [3Ci(/i) + Cal/x) + SCsifx) + din) + SC^ifJ.) + Cein) 
-i^(0, s) {CM + 3C4(/i)) - \g (1, s) {AC, + 4^ + 3C5 + C,) 
-\g (0, s) (C3 + 3C4) + \ (3C3 + C4 + 3C5 + a) , (9) 

where vnc = rric/mb , s = p'^/ml, and 

2 4 2 

^ (^) = -g^^ - 3^^2 (s) - - In {§) In (1 - s) 

5 + 4s , , ., 2s(l + s)(l-2s) , ,., 5 + 9s-6s2 

7 rln 1-s ^ ^ ^Ins +— — 10 

3(1 + 2s) ^ ^ 3(l-s)'(l + 2s) ^' 3(l-S)(l + 2s) ^ ^ 

represents the O {as) correction from the one gluon exchange in the matrix element of Og, 
while the function g {ttIc, s) arises from one loop contributions of the four-quark operators 
Oi-Oq, whose form is 

,. ., 8, _ , 8 4 

g{mc,s) = -- In {mi) + — + -yi 



--{2 + y,) ^\l-yM 9(1- y.) In V,^ [ - t vr 



\/|l-l/*l 



+e(i/,-l)2arctan^i==l , (11) 

where yi = Amf/p'^. 

The Wilson coefficients Cg receives also long distance contributions, which have their 
origin in the real cc intermediate states, i.e., J/ip, ip', ■ ■ ■. The J ftp family is introduced by 
the Breit-Wigner distribution for the resonances through the replacement ([4-7,32]) 

/ - -\ / - -\ 37r .^-^ my r(Vi —> £"•"£") ,.^, 

g {m^, s) ^ g {m^, s) —k }^ ^^^ '_2 ^ , .•___r__ ' ^^^^ 
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em Vi=J/i>i,i>', 

where the phenomenological parameter k = 2.3 is chosen in order to reproduce correctly 
the experimental value of the branching ratio (see for example \^M ) 

B{B -^ j/i^x -^ xtr) = B{B -^ j/4jX) b{j/4j -^ xtr) . 



The effective short-distance Hamiltonian for b -^ si~^i decay [28-31] leads to the QCD 
corrected matrix element (when the s quark mass is neglected) 
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where p"^ is the invariant dilepton mass. 

After obtaining the matrix element for b -^ si^£~ transition, our next task is, starting 
from this matrix element, to calculate the matrix element of the A;, -^ Ai~^i~ decay. It 
follows from the matrix element of the b — >■ si~^i~ that, the matrix elements (A|s7^(l — 
75)6|A;,) and {A\sia^^py{l + 'y^)b\Ah) have to be calculated in order in order to be able to 
calculate the matrix element of the exclusive A;, — > A£+£~ decay. A lot of form factors are 
required for a description of this decay. However when the heavy quark effective theory 
(HQET) has been used, the heavy quark symmetry reduces the number of independent 
form factors for the baryonic transition Ag —>■ light spin-1/2 baryon, only to two [Fi and 
F2), irrelevant to the Dirac structure of the relevant operators (for more details see p3| ) 



(A(p, s) \srb\ Abiv, s')) = u^ip, s){Fi{pv)+ j^F2{pv)}ruA^{v, s') 



(14) 



where v is the four-velocity of A^, F is an arbitrary Dirac structure (in our case F = 7^(1—75) 
and ia^up^i^l + '-y^)). The form factors Fi and F2 for the A;, — > Ai^i~ decay are calculated in 
framework of the QCD sum rules approach in [34| . So the matrix element of the A^ -^ Ai^i~ 
decay takes the following form: 
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Using Eq. (15) and summing over polarization of the final leptons and averaging over 
polarization of the initial A^, we get the following result for the double differential decay 
rate (the masses of the final leptons are neglected and all calculations are performed in the 
rest frame of the A^ baryon) 
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where sa is the spin vector and n is the unit vector along the momentum of the A baryon, 
z = cos6 and the functions A{t) and B{t) are expressed as 
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where r = m/^/mx^ and t = E/rrix^, respectively. It should be noted here that A{t) and 
B(t) were calculated in SM in [Q but our results do not coincide with theirs, especially 
on B{t). Integrating Eq. (16) over t, the differential decay rate can be rewritten in terms 
of the polarization variable a as 

^ = ^{l + {sl-n)a}, (18) 
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and a is the asymmetry parameter, whose form is given as 
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where the integration limits are determined by 
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r < t < - I 1 + r' 
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As we noted previously, in model III a new phase appears. Therefore we would expect 
larger CP violation compared to the SM model prediction. The CP violating asymmetry 
between A^ -^ Ai^i^ and A^ -^ Ai^i~ decays is defined as 

dT df 

dt dt 

The differential widths of the Aj, —^ Kt^l~ and A^, -^ M^l~ decays can easily be obtained 
from Eq. (16) by integrating over z. Hence the CP violating asymmetry takes the following 
form: 

AcP = -^-^^^lraC^/nraa^"{F^it - r') + 2r(l - t)F,F, + F', [{t - r') - 2t(l - t)] } ,(22) 

where 

In derivation of Aqp the following representation of Cg and Cj have been used 

Cl^^ = ReCg + ilmCg , 

C^ff = ReCr + ilmCj , (23) 



and following [^ we assume that the form factors are real. It should be noted that since 
ImCy = in models I, II, and SM, the CP asymmetry is zero (or suppressed very strongly), 
which is one essential difference among the model III and models I, II and SM. 

3 Numerical analysis 

In the present work we have considered three different versions, namely models I, II and III 
of the 2HDM. For the free parameters Xi,b and Xu of model III, we have used the restrictions 
coming from B -^ Xs'~f decay, B^-B^ mixing, p parameter and neutron electric-dipole 



moment [2(;], that yields \\bb\ = 50, |Att| < 0.03. Similar analysis restricts the value of 



tan/3, which is the free parameter of model I and model II, to ||3^, |3^ 

0.7 < tan/3 < 0.6 (m^±/l GeV) , 
(the lower limit of the charged Higgs boson mass is obtained to be mfj± > 200 GeV in 



The values of the main input parameters, which appear in the expressions for the branch- 
ing ratio and Aqp are: m^ = 4.8 GeV, m/^^ = 5.64 GeV, m^ = 1.116 GeV, rric = 1.4 GeV. 
The values of the Wilson coefficients are, Ci = -0.249, C2 = 1.108, C3 = 1.112 x lO'^, 
C4 = -2.569 X 10-2, C5 = 7.4 x 10-^ Cg = -3.144 x lO'^. As has already been noted for 
the form factors that are needed in the present numerical analysis, we have used the results 



of the work ||33 . 





Decay width V x 10^^ (in GeV) 


mH± {GeV) 


Model I 


Model II 


Model III 


100 


6.12 


6.13 


6.09 


250 


6.10 


6.10 


6.08 


400 


6.10 


6.10 


6.08 


1000 


6.09 


6.09 


6.08 



Table 1: 



In Fig. (1) we present the dependence of the differential width of the A;, -^ Kt^t~ decay 
on t at tan/9 = 1.5 and at mu± = 250 GeV for the models I and II. In this figure we also 
depict the dependence of the same differential decay width at lAiiAbf,! = 1.5 and at the value 
of the phase angle = for the same value of the charged Higgs boson mass. In both cases 
the long distance effects are taken into account. 

The values of the decay width in three different models of the 2HDM for different 
choices of the values of the charged Higgs boson mass is listed in Table 1. It follows 
from this table that in all three models the charged Higgs boson contribution to the decay 
width is negligibly small for the values of the \u and \bb (or tan j3) , which lies within the 
experimental bounds. Moreover it should be noted that if the long distance contributions 
{J /ip resonances) are neglected, the decay width becomes two order of magnitude smaller, 
i.e., T.hortdist. ~ 8 X 10~^9 GeV. 

The dependence of the asymmetry parameter a on the charged Higgs boson mass mH± 
and phase angle in model HI is presented in Fig. (2) when long distance effects are taken 
into account. We observe from this figure that the asymmetry parameter Acp increases 
in magnitude as the mass of the m//± increases. This is due to the fact that the decay 
width increases as mH± increases. It is also observed when for < < tt as modulo a 
increases and decreases for tt < < 2ti. This can be explained by the fact that in the region 
O<0<7r(7r<0< 27r), the charged Higgs boson contribution to the SM is constructive 
(destructive). 

For a comparison we present the asymmetry parameter a in SM and 2HDM with and 
without long distance contributions, at mH^ = 250 GeV and = 0. 



a 



a 



short 
2HDM 
long 
2HDM 



-0.48 , 
-0.54 , 



«Im* = -0.50, 
«Im* = -0.54. 

As we have noted earlier, in model III a new phase appears in Xu^bb vertex which is 
embedded in Cj term. As a result interference of the imaginary parts of Cj and Cg 
can induce CP violating asymmetry. 

In Fig. (3) we present the dependence of CP asymmetry on t and the phase angle in 
model III. It is observed from this figure that in the resonance region \Acp\ — 4% , and 
far from resonance region \Acp\ — 1.5%. This is a very useful information in establishing 
model III, since in models I, II and SM Acp is practically zero due to the fact that in all 
these models the Wilson coefficient Cj is real. 

In conclusion, we investigate the A;, -^ Ai~^i~ decay in general 2HDM, in which a 
new extra phase is present. It is shown that investigation of the CP asymmetry which is 
attributed to the differential decay width differences, can give unambiguous information 
about model III, since in this version CP asymmetry can be quite measurable, while at the 
same time CP asymmetry in models I and II are highly suppressed. 
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Figure captions 

Fig. 1 The dependence of the differential width of the A^ -^ Ai^i' decay on t, for three 
different versions of the 2HDM, at rriHi = 250 GeV. The free parameter of models I and 
II is taken tan/9 = 1.5 and for model III we choose \\bb^tt\ = 1-5, |At(| = 0.03. Dotted line 
represents model I, dash-dotted line represents model II and solid line represents model 
III, respectively. 

Fig. 2 The dependence of the asymmetry parameter a on m//± and the phase angle 
for model III. 

Fig. 3 The dependence of CP asymmetry parameter on the dimensionless parameter t 
and the phase angle at m//± = 250 GeV, for model III. 
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